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Abstract 

We derive a number of solution for one-dimensional dynamics of relativistic magnetized plasma 
that can be used as benchmark estimates in relativistic hydrodynamic and magnetohydrodynamic 
numerical codes. 

First, we analyze the properties of simple waves of fast modes propagating orthogonally to the 
magnetic field in relativistically hot plasma. The magnetic and kinetic pressures obey different 
equations of state, so that the system behaves as a mixture of gases with different polytropic 
indices. We find the self-similar solutions for the expansion of hot strongly magnetized plasma into 
vacuum. 

Second, we derive linear hodograph and Darboux equations for the relativistic Khalatnikov po- 
tential, which describe arbitrary one-dimensional isentropic relativistic motion of cold magnetized 
plasma and find their general and particular solutions. The obtained hodograph and Darboux 
equations are very powerful: system of highly non-linear, relativistic, time dependent equations 
describing arbitrary (not necessarily self-similar) dynamics of highly magnetized plasma reduces 
to a single linear differential equation. 



I. INTRODUCTION 



Expansion of plasma into vacuum is a basic problem in fluid mechanics that has a 
wide range of applications from heavy nuclei collisions to astrophysics. In nuclear physics, 
Belenkij and Landau \X\ used the hydrodynamical approach to study multiparticle pro- 
duction in heavy ion collisions. A head-on collision of two highly relativistic nuclei cre- 
ates a relativistically-compressed hot layer of quarkgluon plasma that expands quasi-one- 
dimensionally [2]. On a very different scale, a wide variety of astrophysical objects like jets 
from Active Galactic Nuclei (AGN) [3j, Gamma Ray Burst (GRB) jl] and pulsar winds and 
magnetospheres of a special type of neutron stars - magnetars [5]- may contain relativistic 
strongly magnetized plasma, in which the energy density of the magnetic field dominates over 
the matter energy density (including kinetic and thermal energies), B 2 > pc 2 , P. During 
explosions the strongly magnetized plasma created by the central source suddenly expands 
into a surrounding low density medium. In case of magnetar flares, the initial dissipation 
event creates relativistically hot, strongly magnetized fireball, somewhat analogous to Solar 
coronal mass ejections [6]; the fireball accelerates to relativistic velocities [?]. In long GRBs, 
when a hot magnetically dominated jet reaches the surface of the star it breaks into low 
density medium [U [8]. Similar dynamics may occur in non- stationary outflows in AGNs as 
well [9]. 

In all the above mentioned cases, it is expected that at some distances from the source, 
the magnetic field is dominated by the toroidal component, while motion is preferentially 
radial. Thus, velocity is nearly perpendicular to magnetic field - this type of motion is 
sometimes called transverse magnetohydrodynamics. Equations of transverse MHD reduce 
to fluid equations, but with a complicated equation of state [10] . Qualitatively, such plasma 
behaves like a mixture of fluids with different adiabatic indices, T = 2 for the magnetic field, 
and some T for the kinetic pressure, usually taken to be T = 4/3. Thus, many results of 
fluid dynamics, which assume a single adiabatic index, become invalid. 

Numerical investigation of these phenomena requires the use of relativistic MHD codes 
that can handle high magnetization and high kinetic pressures exceeding the rest mass 
density. Exact, explicit non-linear solutions of fluid equations, and especially relativistic 
MHD equations, are rare. Yet they are important for benchmark estimates of the overall 
dynamical behavior in numerical simulations of relativistic flows and strongly magnetized 
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outflows in particular. In ^TT] we find analytical expressions for self-similar expansion into 
vacuum of a hot magnetized plasma, considering both the Newtonian and relativistic cases 
with arbitrary ratios of kinetic and magnetic pressures to rest mass density. 

In addition, at later times, when the whole initial state of plasma is affected by the 
expansion, the expansion dynamics becomes non-self-similar. A classical related problem is 
then the expansion of a slab of finite length. 

Arbitrary isentropic one-dimensional motion of a fluid is fully integrable. Mathematically 
this is achieved by exchanging the independent variables {t, x} and dependent variables 
{p,v}. As a result a system of two non-linear equations (of mass conservation and Euler's 
law plus assumed isentropic equation of state) is reduced to a single linear equation for the 
(Legendre-transformed) Bernoulli (Khalatnikov) potential [e.g. [TT]. This is always possible 
if the coefficients in these equations do not depend explicitly on time and coordinate. In 
the case of isentropic fluid, the equation of state can be correspondingly inverted. This 



transformation is called the hodograph transformation [llj. In £IV we perform a hodo- 
graph transformation for relativistic cold magnetized plasma and derive the corresponding 
hodograph and Darboux equations. 



II. RELATIVISTIC EXPANSION OF HOT MAGNETIZED PLASMA INTO VAC- 
UUM 

Consider a relativistic one-dimensional flow of hot magnetized fluid along x direction, 
carrying magnetic field in perpendicular direction (so-called transverse MHD) The governing 
of the equations are [TT] 

dt(-Yp) + d x (yPp) = 
d t T m + d x T 0x = 
dtTo x + d x T xx = 
T 00 = l 2 w - P - B 2 /2 
T 0x = l 2 (ivo 

T xx = 1 2 (3 2 w + P + B 2 /2 

£ = p + 3P + B 2 /2, w = p + 4P + B 2 (1) 
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where B is the proper (plasma frame) magnetic field divided by yAn, {3 and 7 are fluid's 
velocities in terms of the speed of light and Lorentz factors, S is the proper energy density, p 
is proper mass density and w is proper enthalpy. The adiabatic index for the kinetic pressure 
is assumed to be T = 4/3. 

Let us consider a semi-space x < occupied by homogeneous hot plasma with density 
Po, magnetic field B and pressure P . The initial state is assumed to be homogeneous. (At 
time t = a barrier located at x = is removed; plasma starts expanding in the positive 
x direction, while a rarefaction wave is launched in the negative x direction. We introduce 
two dimensionless parameters describing the initial plasma magnetization and the ratio of 
kinetic to magnetic pressures: 

a = BI/pq 

B = 2P /B 2 (2) 

Parameter a [12] measures the importance of relativistic effects of magnetic fields: for a > 1 
the Alfven velocity approaches the speed of light. Parameter B is the conventional plasma 
beta parameter: the ratio of kinetic and magnetic pressures. 

The relevant speed of propagation of disturbances is the fast magnetosonic speed Cf 

3B 2 + AP 2 2 _ (3 + 2B)a 



° 2f 3(£ 2 + 4P + p) V ' A + ° s 3( 1 - a( I - 28) ) 

V A = 



.2 a 



1 + ct(1 + 2£) 

s 3(1 + a + Ba) {) 

where va is Alfven velocity and c s is sound speed. Note that for a = 1/2, the fast speed is 
independent of B and equals Cf = l/y/3. In this case the contribution of kinetic pressure to 
total pressure is compensated by its contribution to effective mass density. 

For transverse MHD (when the motion is perpendicular to the magnetic field), the in- 
duction equation and matter conservation require B/Bq = p/po- Assuming that all the 
quantities depend on the self-similar variable i] = x/t, Eqns. ([T]) give 

7 2 (r/ -2(3 + V (3 2 ) (4P + Pl (p + B 2 oPl )) f3' + 

(j]/3(p + 2B 2 2 Pl ) - (3 2 ( Po + B 2 Pl ) - B 2 Pl ) p[ - (1 - A V (3 - 3f3 2 )P' = (5) 



where p\ = p/po and primes denote derivative with respect to rj (note, that if the initial 
state is not homogeneous, the self-similar solutions may depend on combination x a /t). 

Eliminating pressure and magnetic field in the initial state in favor of a and B, P — 
P (p/po) 4/3 , = (Ba/2)p , B = ^/apo, we find 
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1 — 77' 

— 2T a Pi 
T 



1 + 4/5/7 - 3 V 2 + p 2 {rf - 3) )Bap\ 



1/3 







(6) 



Changing to Doppler factors 



'1 + /3 
1-/3 

! l + r] 
1 — 77' 



Eqns. (|4p| become 



0? - - tffiopi + 1$- m + 6*) Bap 1 '" = 



(7) 

(8) 
(9) 



It is more convenient to change to a new variable U = ^p\o and redefine parameter 
B\ = (2/3)/3cr 2 / 3 (in this case the resulting equations (10 11 ) depend only on one parameter 
B\\ also, in the cold limit U becomes a four- velocity of Alfven waves U = — (3\ [II]). 

In terms of variables (3 — U, Eqns (l8p|) become 



(10) 
(11) 



05 



dU 



Previously, Lyutikov [T3] found simple analytical solution of Eqns. (10 11) for simple 



waves in cold magnetized plasma expanding into vacuum. If initially the plasma is at rest, 
and the Alfven velocity in the unperturbed medium va,o is given by corresponding Doppler 
factor 5a,o = y (1 + iU,o)/(l — v A,o), m the limit B\ = 0, Eqns. (10-12) have solutions 



(12) 



In case of hot plasma, B\ 7^ 0, Eq. (10) can be resolved for S 



S 2 1 + AB.U 2 / 3 + 2U 2 ± 2U 1 / 3 J3B 2 1 U 2 / 3 + U^ 3 + U 1Q / 3 + B 1 (1 + AU 2 ) 
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2BXU 2 / 3 



f(U). (13) 



This gives a general solution for the Doppler factor in terms of initial parameters B\ 
(2/3)i3a 2 / 3 and local density U = . 



Eq. (11) then becomes 

dU f(i + f 2 W 



d\n5 v 2f(l-P)-(l + P)Uduf 



or, using the explicit form of /, Eq. (13), 



dlnSr, 12B 2 U 2 ' 3 + W 4 ' 3 + Bi(7 + 16U 2 



(14) 



(15) 



du 3(1 + 2B 1 U 2 / 3 U 2 / 3 ^f(B 1 + t/ 4 /3)(i + 3B 1 U 2 / 3 + U 2 ) 
Changing independent variables 5 V — >■ U, we find: 

US, = - 1° 1^ + 9^ + 8,(7 +16^) 

7 3(l + 2B l2 / 2 / 3 )y 2 /3 v /(B 1 +y 4 /3)(l + 3B iy 2/3 + y 2) 



Eq. (16) gives a general solution for simple fast waves in relativistic magnetized fluid with 



arbitrary ratios of magnetic and kinetic pressure to rest mass, with kinetic part of the 



pressure obeying adiabatic law with T = 4/3. Eq. (16) expresses implicitly the density 
U = y/plcr in terms of the self-similar variable rj = x/t. 



The lower limit of integration in Eq. ( 16 ) can be found from the conditions on the front 
characteristics propagating into undisturbed plasma. The rarefaction wave propagates into 
the undisturbed medium with the local fast velocity (see Eq. (pi)): 

2 B 1 a 1 / 3 + a 

Cf '° ~ 1 + 3BktV3 + a {U) 

This corresponds to 



"770 



1 + c /,o 



^•=5(^-4.) (") 

(note the signs of velocities in the definition of 5 m : the rarefaction wave propagates in the 
direction opposite to the flow). At the front of the rarefaction wave 5 V = 5 Vo . Thus, a 
particular solution corresponding to expansion into vacuum of a medium initially at rest is 



6^ ru 3(1 + 2ff 1 y 2 /y/ 3 v /(ff 1 + yV3 )(1 + 3B lV 2 / 3 + y 2 ) 
n 5 V0 ~ Juo \2B 2 y 2 / 3 + V/3 + B x (7 + 16y 2 ) V { } 



Eq. (19) gives a solution for relativistic expansion of hot magnetized plasma into vacuum, 



see Fig. [TJ Corresponding velocities are plotted in Fig. [2} 



One can verify that Eq. (19) reproduces the known result for cold magnetized plasma 
For zero kinetic pressure, B\ = 0, we find 

f(U) = 5? 

Su = Of 73 (20) 

in agreement with [1] 



For zero kinetic pressure, B\ — > 0, relations (|19||20|) imply that the front of the rarefaction 
wave is located at tjrw = 



'cr/(l + cr), while for large kinetic pressure, B\ — > oo, the front 
of the rarefaction wave is located at riRw = — 1/\/3. The vacuum interface in the limit 
£>i — > oo approaches the speed of light. 
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FIG. 1. Four- velocity of the relativistic self-similar expansion for different values of B\ = 
(2/3)£ct 2 / 3 , B = 2P /B$ . Plotted is the value of U = yjpio as a function of the self-similar 
parameter for various values of the kinetic pressure parameter B\ and magnetization. Left Panel: 
a = 1, Right Panel a = 1/2. The upper curve, corresponding to B\ = 0, is, in fact two coincident 
curves, plotted using the known explicit solution for relativistic expansion of magnetized cold gas 



into vacuum, Eq. (12), see also [T3], and the implicit solutions (16) for zero kinetic pressure. In 
this case the front of the rarefaction wave is located at rj = —\Joj(\ + a). For a = 1/2, the front 
of the rarefaction wave is located at r\ = — ^Ja j{\ + a) = —l/y/3 

For a fixed values of a > 1/2, the increase of kinetic pressure leads to larger values of 
negative r/nw (smaller absolute values of r/nw) due to effective increase in plasma inertia. 
For a = 1/2, the front of the rarefaction wave is located at fixed i] = — 1/a/3 for any B\. 
Finally, we note that the Riemann invariants in a transverse relativistic MHD, does not have 
a representation in simple functions. 
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FIG. 2. Velocities for relativistic expansion of hot plasma into vacuum a = 1 (Left Panel), a = 1/2 
(Right Panel). 

III. SELF-SIMILAR EXPANSION OF UNMAGNETIZED RELATIVISTIC FLUID 

In this section we simplify the above relations for the case of simple waves in relativistic 
unmagnetized fluid. Adopting a polytropic EoS, we will calculate simple waves in plasma 
with adiabatic index of 4/3, yet for a finite ratio P/p, without assuming that the speed 
of sound equals l/v3. This is an important step since most codes use a single value of 
adiabatic index, but for finite ratios Pj p. 

For a finite ratio of P/p, the speed of sound equals 

° s ~ 3(4P + Po )' (21) 

while Eqns. do]) become 



(V ~ Pf - ^ (l " 3r/ 2 + 477/3 - (3 - r? 2 )/? 2 ) p{ /3 = (22) 
6p ^ / 

P\ = j^f 7 2 Px/3' (23) 

In case of negligible density p 0, this immediately gives the ultra- relativistic limit 
= 1/V3, 

p_ 2r ] ±V3(l-r ] 2 ) 
3 — i] 2 

(24) 



2- V3 



The front of the rarefaction wave is located at i]^/ = — l/v3, while on the vacuum side 
expansion proceeds with the speed of light. 

In the ultra-relativistic limit the flow lines dx/dt = (3 are given by 

V3/2 



t 

to 



1 



l+T] 



y/1 — T] 2 \ 1 — T] J 

The characteristics satisfy dx/dt = ((3 + c s )/(l + (3c s ), which gives 

1/V3 



(25) 



t 

to 



(the other characteristics are straight lines). 



1 +T] 

l-r l/ 



(26) 



In case of finite density, introducing Doppler factors, Eqns. (22 23) give 

1/3 



3 l/3 



Pi 



0? - * 



pi 



Po 



2 \^-6l-6^ n P 



(27) 
(28) 



Which can be integrated 



(3 + 2VS)C - (3 - 2V3)5l V * 



(29) 



The constant of integration C in Eq. (|29j) can be found from the condition that at the front 
of the rarefaction wave, propagating with velocity — c Sj0 and located at i] = — c Si0 , <5»?,o = 
+ r] )/ (1 — 770) the fluid is at rest, /3 = 0, 8p = 1. This gives 



C 



^+(3-2^)5.; 



T/.O 



\/3 - (3 + 2^3)^ 



(30) 



T/.O 



Equations (27, 29, 30) give an analytical solution to the problem of self-similar expansion of 
fluid with T = 4/3 into vacuum, valid for arbitrary ratios of kinetic pressure to mass density, 
see Fig. |3j 



IV. NON-SIMILAR EXPANSION OF COLD MAGNETIZED PLASMA: RELA- 
TIVISTIC HODOGRAPH AND DARBOUX EQUATIONS 

Relativistic hodograph transformation was derived by Belenkij & Landau pQ for poly- 
tropic fluid. Here we first re-derive the corresponding equation for cold magnetized plasma 
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FIG. 3. Self-similar expansion of relativistic unmagnetized fluid with adiabatic index T = 4/3. 



Top to bottom Po/po = 0.1, 1, 10, 100. The curve Po/po = oo is given by Eq. (24). The vacuum 
interface propagates with the speed of light. In the ultra-relativistic limit, Po/po —> oo the front 
of the rarefaction wave is located at t\rw = — l/v3- 

and then transform it to a normal form, where Riemann invariants are taken as independent 
variables. In case of regular fluids the internal energy and the corresponding hodograph 
equations are denned in term of temperature, which is zero in the case of cold magnetized 
fluid. As we will see below, the role of temperature is taken by the proper enthalpy. 

The relativistic hodograph transformation is achieved by introducing the Khalatnikov 
potential cj) [TT] 

j/3w = d x <p 
jw = -d t (f> 

w = ^ (31) 

In the non-relativistic limit the corresponding equations are the condition on potential flow 
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and the Bernoulli equation. Below we drop the tilde sign over the enthalpy: w is then the 
enthalpy per unit mass: 

The differential of the Khalatnikov potential is 

dip = d x <pdx + d t (f>dt = -yfiwdx — jwdt (33) 

Next we employ Legendre transform of <fi with respect to variables {x,t}. The transformed 
potential x becomes 

X = <t> — j(3wx + jwt 

d x = (t - Px)idw -{x- /3t)wy 3 d/3 (34) 

Variables t and x are then expressed from the potential x as 

, dx _ J^dx 
dw •yw d(3 

x = pt-^% (35) 
j 6 w op 

These equations can be modified if we introduce rapidity /3 = tanh r: 

dx sinh r dx 



t = cosh r 



dw w dr 



. dx coshrdx 

x = smhr— — Mo) 

dw w dr 

Comparing with Ref. [1] Eq. (4.12), we see that in cold magnetized plasma proper enthalpy 
w plays a role of a temperature, while the speed of sound c s - Alfven velocity - is given by 

c -= v *=%=wh =1 - 1/w (37) 

The relativistic hodograph equation is then obtained from the continuity equation by 
transforming to independent variables r and w. 

d 2 rX ~ wd wX + (1 - w)wd 2 wX = (38) 
This is relativistic hodgoraph equation for one-dimensional motion of cold magnetized 



plasma. Eq. (38) reduces to the one obtained by Belenkij and Landau for hot fluid with a 
substitution w = (p + B 2 )/ p — > T. A general separable solution of Eq. (38) is 



x = e Cr w 2 F 1 (l-C,l + C,2,w) (39) 
11 



where, we remind, r is the rapidity, = tanhr, and 2-^1 is the hypergeometric function. 



Let us transform the hodgoraph equation (38) taking the Riemann invariants as indepen- 
dent variables. The Riemann invariants are 



Ji 
J2 



log 8\8p 



log 



S 2 

°A 



(40) 



where 



1 1 + Pa 
I -Pa 

+ 



1 + yjl - 1/w 



l/w 



1 + tanh r 
1 — tanh r 



(41) 



are corresponding Doppler factors. 



Using Eqns (41), the Riemann invariants become 



Ji = 2arctanlW 1 — l/w + r 



J2 = 2arctanhy 1 — l/w — r 

Thus, the rapidity r and the proper enthalpy are 

J\ — J2 



(42) 



w = cosh 



2 J\ + J2 



(43) 



Changing the independent variables in the hodograph equation (38), we derive the rel 



ativistic Darboux equation for the Khalatnikov potential as a function of the Riemann 
invariants: 

1 /), V4- flr V 

(44) 



AX +7 ^**^ =0 
1 2A 4 sinh^a 



In the non-relativistic limit J\ -> (3 + 2(3 A , J 2 -> 2^-/9, sinh (Ji + J 2 )/2 (Ji + J 2 )/2 and 



Eq. (44) reduces to the non-relativistic Darboux equation for a fluid with adiabatic index 



of T = 2, Eq. (B3) 



Using constancy of the Riemann invariants on the characteristics, one finds (c/., Eq, (45)) 



dj 2 x = tanh ( 1 - — - ) OjJ 
d j2 x = tanh (^^^ ) 0.,J 



(45) 
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System (45) can be written as a single equation for time 



3 sech 2 d Jx t + sech 2 (^f^) dj 2 t 



Ji-3J 2 



(46) 



4 tanh (^f^) - tanh (- 
This is the relativistic analogue of the Darboux equation for the time variable. Solutions 



of Eqs. (46 45) give time and spacial coordinate as functions of two Riemann invariants. 
These solutions can then be inverted for Ji^(x,t). 



V. ID EXPANSION OF A SLAB OF GAS: SOLUTIONS OF THE HODOGRAPH 
EQUATION 

Consider a slab of gas initially at rest, occupying region < x < L and expanding 
into vacuum x > L. At x = there is an impenetrable wall, Fig. |4j In the initial state 
the Riemann invariants are J^o = J2,o = 2c Sj o where c Sj o = ^a,o is the sound (Alfven) 
velocity in the initial state. The boundary conditions for this problem are [Hi problem after 
parag. 105] : zero velocity at the wall and constancy of the first Riemann invariant on the 
characteristics that leaves the wall at the moment of reflection: 



dx 

d(3 



= 

x(Ji = Jifl) = (47) 



Numerical solutions for Alfven and sound speeds are found and plotted in both the 
non-relativistic (as discussed in Appendix |B| and relativistic case. In both cases find- 
ing a numerical solution involves inverting a pair of functions of the Riemann invariants, 
{t ( Ji, Jq) , x (Ji, J2)}, to obtain functions of position and time giving the Riemann invari- 
ants, {Ji(x, t), J2(x, t)}. In the non-relativistic case the initial functions {t ( Ji, J2) , x ( Ji, J2)} 



are functions represented by Eq.'s (B8) and (B9) whereas in the relativistic case these func 



tions are determined using numerical methods. The relevant quantities, Alfven and sound 
speed, can then be expressed as functions of position and time. The numerical solutions 
in both relativistic and non-relativistic cases involve separate calculations for the regions 
of expansion affected by the reflection of the initial rarefaction wave at the position of the 
wall, x = 0, and the regions unaffected by this reflection (as well as the boundary between 
these regions). The former region will be referred to as the "complex" region and latter as 
the "simple" region. 
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B. Numerical Solutions of Relativistic Expansion 



Initial conditions in the relativistic case are given, by initial Alfven speed (3a(x, 0) = 
tanh(l), initial velocity (3(x,0) = 0, the plasma occupies the region < x < tanh(l) with 
an impenetrable wall at x = and the plasma expands into vacuum at t — 0. The initial 
conditions are such that reflection of the initial rarefaction wave occurs at t = 1. The Initial 
Alfven speed, tanh(l) ~ .762, was selected as a speed somewhat near the speed of light that 
gives convenient initial values for the Riemann invariants: 2 = J\ o = </2,o- 



In the case of relativistic expansion, no simple analogue to Eq. (B8) that satisfies Eq 



(46) and gives time as a function of the Riemann invariants in the complex wave region is 
readily available. Therefore, in the complex region, time must be numerically calculated 
as a function of the two Riemann invariants, Ji j2 . The first step in this calculation is 
accomplished by numerically solving for a function satisfying the adjoint differential operator 
corresponding to Eq. (73) subject to boundary conditions as outlined in [151 Ch. 5, Sec. 2, 
(Eq.'s 3 and 3')], obtaining a solution R [J 1( J 2 ; a, (3}. In the present calculations this function 
is solved numerically after specifying boundary values based on a and /3. Exploiting the 
symmetry of the Riemann's function, R [J ly J 2 ; a, f3] [T5J Ch. 5, Sec. 2, (Eq.'s 3 and 3')],, 



one obtains an expression for time given as a function of J\ and J 2 satisfying (46) : 



t [Ji, J 2 ] = R ^tanh" 1 (/3 Aj0 ) , 2tanh" 1 (/3 Aj0 ) ; h, J 2 ] (48) 

The initial value of both Riemann invariants, 2 tanh -1 ((3a,o), corresponds to the initial 
conditions of /3 t=tQ = and 0A,t=t o =Pa,o- Position can then be calculated as a function of 



the Riemann invariants by integrating by parts Eq. (45). 

As a next step towards a numerical solution, J 2 is calculated as a function of J\ and time 
by identifying the values of J 2 corresponding to a specified time through an interpolation 
method similar to the one used in the non-relativistic case, resulting in the calculation of 
a function J 2 (Ji,t). In order to determine the numerical solutions in the complex region 
for a given time, a collection of points of the form ( J 1; J 2 (J x ; t) , x [J%, J 2 ( J\\ £)]), where we 
have emphasized only a parametric dependence on time, are calculated based on a number 
of J\ values in the range (Ji, x =o(t), 2). This range represents the J\ values in the complex 
region and Ji, x =o(t) represents the minimum value encountered at the point x = subject 



to the condition t = R 



2 tanh" 1 (p A>0 ) , 2 tanh" 1 (/3 Aj0 ) ; Ji, x=0 (t), Ji, x=0 {t) 
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. This condition 



Jl,x=Q — J2,x=0 
2 



is based on the requirement that /3(0, t) = = tanh ( 
Having obtained a sufficient collection of points, the solutions 

Ji(x,t) - J 2 {x,t) 



and hence J\ 



2 = 



J. 



2,x=0- 



f3(x,t) = tanh 
/#a(:e, £) = tanh 



Ji(x,t) + J 2 (x,t) 
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(49) 



can be numerically calculated by interpolation. The functions 0(x,t) and 0A(x,t) in the 
simple region of the expansion are known analytical functions of x and t, more specifically 
of a single variable rj = x =^^ given in Ref. [H] as 

1 



~ 1 + 5 A0 4/3 V /3 
0a(v) = 1-2 



A 2/3 



^ 2 / 3 + 5ao 4/3 



(50) 



The boundary, x;,(t), between the simple and complex regions of expansion is determined as 
a function of time in terms of r/b = Xb ~ x ° by the equation 



<Uo 4 



l+5n 2 



'Vb 



(51) 



Combining the interpolated complex region solutions and the simple region solutions given 
in Eq.'s (50), complete solutions are plotted in Figure Q. 

Previously, the problem of non-self-similar expansion of magnetized plasma was consid- 
ered numerically in Ref. [I]o\. Contrary to the initial claim in Ref. [15], the presence of the 
wall is detrimental to acceleration, as can be seen from Fig. |4j Magnetic pressure-driven 
acceleration proceeds most efficiently while a given fluid element is causally disconnected 
from the wall, during the self-similar stage discussed in Ref. [1] 



VI. DISCUSSION 



In this paper we derived a number of analytical results for one- dimensional expansion of 
magnetized gas into plasma. First, we found the self-similar expansion into vacuum of a 
hot magnetized plasma. In this case the total pressure has contribution both form magnetic 
field and from kinetic motion of particles. These two contribution obey different equations 
of state, so effectively, we considered relativistic self-similar expansion of a mixture of gasses 
with different adiabatic indices. 
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FIG. 4. Non-self-similar one-dimensional relativistic expansion of a slab of magnetized plasma, 
initially occupying < x < 1 and limited by a wall at x = 0. Expansion proceeds in the positive 
direction. Plots of (a) velocity and (b) Alfven speed, along with related 7-factors (c-d) with respect 
to position measured from the wall at x = for times t =1, 2, 5, and 10. Velocities are measured 
as fractional multiples of c and position, x, is measured in the same units as time multiplied by the 
speed of light. Initially the plasma occupies the region < x < Paq (Note: Pa = tanh(l) in the 
present calculations) and expands into the vacuum occupying x > (3ao while an impenetrable wall 
is present at x = (resulting in the reflection of the initial rarefaction wave at t = 1). The effect of 
finite distance from the wall occurs for times after t = 1 and can be observed as the discontinuities 
of the first derivatives in the velocities and Alfven speeds, marking the separation between complex 
and simple regions of expansion. 

Second, we derived relativistic hodograph and Darboux equations that describe arbitrary 
one-dimensional motion of magnetized plasma perpendicular to magnetic field. The obtained 
resulting hodograph and Darboux equations are very powerful: we reduced a system of highly 
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non-linear, relativistic, time dependent equations describing arbitrary (not necessarily self- 
similar) dynamics of highly magnetized plasma to a single linear differential equation. Using 
semi-analytical methods we calculated evolution of the flow parameters. 
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Appendix A: Non-Relativistic expansion of hot magnetized plasma into vacuum 

In the non-relativistic limit, the equations of one-dimensional transverse MHD read 

d t p + d x (vp) = 

p(d t v + vd x v) = -d x (B 2 /2 + P) 

d t B + d x (Bp) = 

P = P (p/po) r (Al) 



where P is kinetic pressure, B is magnetic field divided by v47r, p is density is v is plasma 
velocity. The kinetic pressure obeys a polytropic equations of state with index T. As the 
initial condition, we assume that at time t — plasma occupies region x < 0, while at 
x > the medium is a vacuum. Initial homogeneous density, magnetic fields and kinetic 
pressures are po, Bq and Pq correspondingly. At time t = a barrier at x = is removed 
and the plasma starts expanding into vacuum while a rarefaction wave propagates into the 
bulk plasma. In the initial state the Alfven and sound velocities are 

v \o = b o/Po 

cjo = rP /p (A2) 
Let us assume that all the quantities depend on self-similar combination rj = z/t. The 



system (Al) then reduces to 

(77 - v)p[ + v'pi = 

(v-v) 2 = v% Q p 1 + cl p T 1 - 1 (A3) 



where p\ = p/po- Eq. (A3) can be resolved for velocity v(px), 

V=T)± V^oPl + ^oPl" 1 , ( A4 ) 

and an equation for p\\ 

d " Pl = ± ^ AI> + < , )P \'M + r) (A5) 
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Eq. ( A5 ) can be integrated 



V = 0^ — ^ + 01^ - 2 (3 _ r)(r _ 1) — p[ 



r -2 



A,0 (3-r)/2 



i 3 - r 7 - 3r 



2' 2(2 -T)' 2(2 -T)' < 
(A6) 



where 2 -Pi is a hypergeometric function. Signs in Eqns (A5) and (A6) correspond to the 



choice in (A4) 



The constant of integration C\ can be found from the condition that at the front of the 
rarefaction wave, which propagates with fast velocity in the undisturbed medium Vf$ and 



is located at tjrw = — f/,o = ~\/ v ao + c s,o; the plasma density is undisturbed p\ = 1: 



Ci 



r-2 



< „ (1 3-r 

"2-Tl 



7-3r 



( 3_r)(r-i)c Si0 z 1 V2'2(2-r)'2(2-r)' < 

In particular, for adiabatic index T = 5/3 the previous relations simplify 

1 



(A7) 



V 

C s ,0 



V 



M 



[2(1 + M\ fi f 2 - (2 + SM^p 1 / 3 )^ + Ml p\ /3 

n (l + ^l,o) 3/2 ~ (1 + ^Vl /3 ) 3/2 
2 -j^*- 



A,0 



c s,0 



(A8) 



see Fig. [5] The front of the rarefaction wave propagates with the fast speed and is located 



at T]rw — — c S) oyl + M\ Q , while the vacuum interface, corresponding to p\ =0, propagates 
with the velocity 

o (l + Mj, ) 3 / 2 -l 

Uvnr = 2 7^ C Sj0 (A9) 



For unmagnetized plasma Eq. (A8) gives 



M 2 



jj = (3 - 4p} /3 )c Sj0 

^vac = 3c s 0) 



(A10) 



while for cold magnetized plasma 



77 = (2 - 2y/pdv A>Q 

V vac = 2v A n, 



(AH) 



in correspondence with the general solutions in media with adiabatic indices T = 5/3 and 

r = 2 HHE]. 



v A,o 2-r 
Pi 
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FIG. 5. Self-similar non-relativistic expansion of magnetized fluid into vacuum, Eq. (A8), T = 5/3. 
Left Panel: the value of the self-similar parameter in terms of fast velocity, r]/vf t o = rj/ ^Jc? s + v\ 
as a function of density p\ = p/po- Right Panel: velocity in terms of sound velocity, v/c s> o as 
function of density. The front of the rarefaction wave is located at p\ = 1, the vacuum interface 
is located at p\ = 0. Different curves correspond to different values of Ma,o = va,o/cs,o = 0, 1, 10 
(top to bottom). 

More, generally, in case of unmagnetized plasma, va,o = 0, and arbitrary adiabatic index 

r, 

(So that the vacuum interface expands with 2c Si0 /(r — 1).) 
The Riemann invariants [TT] are 

J ± = v±[^ = v±2 C -^i (A13) 

J pv f (?y A0 



where p is total pressure, kinetic plus magnetic. Relations (A13) assume that both sound 
and Alfven velocities are non-zero. 

The Riemann invariants J± are constant along the corresponding characteristics C±: 
dz/dt = v ± Vf. We find then that the C_ characteristics are straight lines x = rjt , while 
C + characteristics are determined by 



dz _ / (1 + M 2 ) 3 / 2 y/l + M^ /3 (2 + MV 1 /3 ) , 

dt ~ 2 W hP 1 Cs '° (A14) 



where px(z/t) should be found from Eq. (A8). This gives a transcendental equation for the 
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characteristics: 
dz\ 3 _ ( 



! \ 

-c s , r / + 16((l + M 2 ) 3 -l/9) <? a A 



1 f n » I i 2 ( 1 + M2 ) 3/2 c r? I 16 3(1 + M2)3 + 1 c 2 n I 61 v / miP(3 + 6M 2 + 4M 4 + M 2 



Appendix B: Non-relativistic hodograph and Darboux equations for cold magne- 
tized plasma 

Let us next discuss non-relativistic one-dimensional motion of cold magnetized plasma. In 
this case the induction and continuity equations imply B/p =constant, so that the pressure, 
P oc B 2 is related to density via adiabatic law P oc p r with T = 2. 

For polytropic index of T = 2, the corresponding hodograph equations becomes [HI §105] 



wd w x - d v x + d v x = 

t = — 

w 

x = vt - d * (Bl) 

v 

where w = c 2 s = v \ is proper fluid enthalpy and c s = Va is the sound (Alfven ) speed. 



A general solutions of the equation (Bl) can be easily obtained. For example, a solution 



separable in w,f3 is \ = e a Plo(2ay/w). The main problem with solving the hodograph 



equation (Bl) (or the corresponding Darboux equation (B3)) in a simple application of 
one-dimensional expansion of gas into vacuum, is in finding solutions that satisfy boundary 
conditions, one given on the characteristics and another at a fixed velocity. 



Equation (Bl) is often transformed taking Riemann invariants as independent variables 

m\ 

Ji = 2c s + v 

J 2 = 2c s -v (B2) 

Ji is constant on the characteristics dx/dt —v + c, while J 2 is constant on the characteristics 
dx/dt = v — c. In the initial state J\ = J 2 = 2c s>0 = Ji,o = ^2,0- 
In terms of { J±, J 2 }, the hodograph equation becomes 

dj^x + Af dj if + dj * x =0,Af= 1/2 (B3) 
Ji + J 2 
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Eq. (B3) is referred to as Darboux equation. Darboux equation played an important role in 
the development of fluid mechanics: analysis of the analogue of the Darboux equation lead 
Riemann to the formulation of the theory of hyperbolic equations. 

For integer values of the coefficient M = (1/2) (3 — T)/(r — 1), Darboux equation can be 
reduced to a one-dimensional wave equation. For the case of interest, r = 2, M = 1/2, this 
is not possible. 

Equivalently, using constancy of the Riemann invariants on the characteristics, one find 

3 Ji — J2 n 

°J2 X = 4 °Jat 

d j2 x = J -^^d j2 t (B4) 
On can then write Darboux equation for time variable: 



The main mathematical difficulty in solving Darboux equation (B3) with boundary con- 



ditions (47) is that one boundary condition is given on the characteristics, while the other 
at a fixed values of velocity. For integer values of Af [U] give the solution; e.g. for Af = 1, 



A general solutions of the Darboux equation (B3) can be expressed in terms of the 



corresponding Riemann function (the analogue of the Green's function) [18, Eq.3.71] 



Wl,0 + J 2,0 J \ { J l ~r J 2) W1,0 "T J 2fi) J 



It turns out that the corresponding Darboux equation for time, Eq. (|B5j) with boundary 
condition t = t , the moment of reflection, when J\ = J 2 = 2 can be solved explicitly [see 
[13 Eq. 82.17] 

t — toB(Ji >0 , J 2) o, Ji, J2) (B7) 

(In passing we note that the equations of the one-dimensional fluid motion with adiabatic 
index T = 2 are equivalent to shallow water equation. The corresponding problem of a dam 
break with a finite lock length has been solve in Ref. [20J.) 
This gives 

t(T n - 64 F (3 3 1 (Ji,o - Ji)(J2, - J2) \ mR] 

{ u 2) - ij 1+ .h^.h + .h: 1 U' 2" (^o+ji)(j 2 ,o+j 2 )J (B8) 
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where 2 -Pi is the hypergeometric function. 



Coordinate x then can then be derived from Eq. (B4): 



x = *J±±A t _ I f h td j 2 (B9 ) 

In the numerical solution of the non-relativistic case, the initial conditions are such that 
the velocity, v(x, 0) = 0, the Alfven speed, ca(x, 0) = 1, the plasma occupies the region 
< x < 1, and an impenetrable wall is present at x — 0. These conditions are exactly 
analogous to those given in Ref. [20] for the case of large Froud number. The front of 
the expansion of plasma propagates as XFront(^) = 2t + 1, corresponding to the foremost J 2 
characteristic, with Riemann invariant value J2 = —2, emanating from x = 1 at t = 0. (We 
refer to characteristics for which ^ = v — cas J2 characteristics and characteristics for which 
n| = v + c as Ji characteristics, the names thus referring to the Riemann invariant which 
is constant on the specified characteristic) The linear relationship governing the expansion 
front position as a function of time is obtained by integration of (Eq. 48) . At time t — 1, 
the initial rarefaction wave reaches the wall at x = and is reflected. Boundary conditions 
are such that plasma velocity is at the wall location x=0 (i.e., J\ — J2 = when x — 0). 
The constancy of the Riemann invariant J2 = 2c — v on the characteristic % = t> — c, implies 
the leading rearward propagating characteristic emanating from x = 1 has value J 2 = 2 on 
the path x — 1 — t. Correspondingly, the J\ characteristic emanating from x — 1 at t — 
has the value 2, satisfying the condition J\ — J2 = when x = 0. All J\ characteristics will 
have a constant value J\ = 2 for £ < 1 in the region {1 — t < x < 2t + 1}, satisfying the 
initial conditions of v(x, 0) = and c^(x,0) = 1. For t > 1, the Ji characteristics begin 
to influence the flow in the region < x < Xb(t). The discontinuities in the first derivatives 
of the plots after t = 1 separate the simple and complex region of the flow. The conditions 



determined by eq.'s B4 in the simple region where J\ = 2 are 



dj 2 x = - (6 - J 2 ) <9 j2 t 

x = 1 + ^(2-3J 2 )i (B10) 

We can numerically solve for the position, Xb(t), of the boundary between the simple and 
complex regions. This point is given by the implicit equation 

^)- 1 + a - (2 2 + ~S)^ (B11) 
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The solution for J 2 (and therefore also velocity and Alfven speed) in the simple region, where 



J\ = 2, can also be obtained from Eq. (BIO), which gives 

1 



J 2 (x,t) 



(B12) 



The functions {t (Ji, J 2 ) , x (Ji, J 2 )} are given, in the complex region, by Eq.'s (B8) and (B9) 



The first step in numerically constructing the functions { Ji(x, t), J 2 (x, t)} in the complex 
region is accomplished by numerically constructing a function J2(Ji,t) by interpolating 
values of t(J\, J 2 ) at numerous values of J 2 over the range [— 2, Ji] (this range covers all 
possible values of J 2 in the complex region). Having achieved computationally constructing 
a function J 2 (Ji,t), it is then possible to calculate position as a function of time and J\\ 
x(Jx,t) = x (Ji, J 2 (Ji, t)). A procedure similar to that used in calculating J 2 (Ji,t) from 
t (Ji, J 2) is then be employed to calculate J%(x, t) and thereby J 2 (x, t) = J 2 {Ji(x, t),t). The 
procedure described allows for the Riemann invariants Ji j2 , and hence v(x,t) and c^(x, t) 



via Eq. (40), to be calculated in the complex region. 
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FIG. 6. Non-relativistic expansion of magnetized plasma limited by the wall at x = 0. Plots of the 
plasma velocity, v(x, t), and Alfven speed, c(x, t), as functions of distance from the wall at x = for 



one-dimensional flow satisfying Eq.'s (Bl), (47). Initial conditions are such that gas is initially at 
rest occupying the region < x < 1 with Alfven speed ca = 1. An impenetrable wall is present at 
x = and the gas expands into vacuum at t = 0. Times shown in this plot are: t = 1, 2, 5, and 10. 
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